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Abstract 

We study the sensitivity of neutral-current neutrino-nucleus scattering to the strange-quark 
content of the axial-vector forni factor of the nucleon. A model-independent formalism for this 
reaction is developed in terms of eight nuclear structure functions. Taking advantage of the insen- 
sitivity of the ratio of proton (v, v'p) to neutron (u, v'n) yields to distortion effects, we compute 
ali structure functions in a relativistic piane wave impulse approximation approach. Further, by 
employing the notion of a bound-state nucleon propagator, closed-form, analytic expressions for 
ali nuclear-structure functions are developed in terms of an accurately calibrated relativistic mean- 
field model. Using a strange-quark contribution to the axial-vector form factor of g s A = —0.19, a 
significant enhancement in the proton-to-neutron yields is observed relative to one with g^ = 0. 
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I. INTRODUCTION 



Neutrino physics has established itself at the forefront of current theoretical and exper- 
imental research in astro, nuclear and particle physics. While neutrino-oscillation experi- 
ments, which are presently evolving from the discovery to the precision phase, will remain 
at the center of most investigations, a variety of other interesting (non-oscillation) physics 
topics may be studied in parallel. A prime example of such a paradigm is the recently 
commissioned MiniBooNE experiment at Fermilab. While MiniBooNE's primary goal is 
to confirm the neutrino-oscillation experiment at the Liquid Scintillator Neutrino Detec- 
tor (LSND) at the Los Alamos National Laboratory 1] this unique facility is ideal for the 
study of supernova neutrinos, neutrino-nucleus scattering, and hadronic structure. In this 
contribution we focus on hadronic structure, in part as a response to the Fermilab Intense 
Neutrino Scattering Experiment (FINeSE) initiative that aims to measure the strange-quark 
contribution to the spin of the proton via neutral-current elastic scattering 0. 

A measurement of the spin asymmetry in deep-inelastic scattering of polarized muons 
on polarized protons by the European Muon Collaboration [3( revealed a disagreement with 
the Ellis-Jaffe sum rule in an approach that assumed that only up and down quarks (and 
antiquarks) contribute to the proton spin. This was one of the first indications that hidden 
flavor in the nucleon may play an important role in the determination of the spin struc- 
ture of the proton. Experiment ally, the spin structure of the proton is also accessible via 
parity-violating electron scattering. Unfortunately, large radiative corrections ja, |y] as well 
as nuclear-structure effects (tJ hinder the extraction of strange-quark information. A comple- 
mentary experimental technique that may be used effectively to study the spin structure of 
the proton is elastic neutrino-proton scattering. The advantages of neutral-current neutrino- 
proton scattering over parity-violating electron scattering and deep-inelastic scattering are 
well documented in the literature @> H[. Two notable examples are: (i) the insensitiv- 
ity of the extraction of the strange-quark contribution to the use of (broken) SU(3)-flavor 
symmetry [Ì[ and (ii) the absence of radiative corrections in neutral-current neutrino scat- 
tering @,E3|. The quark structure of the nucleon may be investigated in a particularly clean 
fashion by evaluating matrix elements of suitable quark-current operators between single nu- 
cleon states 0, H, E3 11 1- This is because quark-current operators may be written in terms 



of the fundamental couplings of the quarks to the Z°-boson, which are fully prescribed in the 
Standard Model. Further, the weak neutral current of the nucleon may be parametrized on 
completely general ground in terms of two vector and one axial form factors (an additional 
induced pseudoscalar form factor is present but its contribution vanishes in the limit of a zero 
neutrino mass). In particular, the axial- vector form factor may be split into a non-strange 
contribution, that may be determined from nuclear beta decay, and a strange contribution 
proportional to the fraction of the nucleon spin carried by the strange quarks [12]. Thus, 
the axial- vector form factor is cruciai to understanding the role played by strange quarks in 
determining the properties of the proton and represents the main focus of this contribution. 

A measurement of neutrino-proton and antineutrino-proton elastic scattering at the 
Brookhaven National Laboratory (BNL) reported a non-zero value for the strange form 
factors of the nucleon l| • However, these results must be treated with caution as the value 
of the axial mass Ma and g s A are strongly correlated [12]. Another point of concern in the 
BNL experiment was that 80% of the events involved the scattering of a neutrino off carbon 
atoms and only 20% where from free protons. Before any finn conclusion may be reached, 
it is therefore necessary to understand nuclear-structure effects. Unfortunately, scattering 
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off a nucleus introduces its own complications. These include: (i) questions concerning 
finite-density effects, such as possible modifications to the nucleon properties in the nuclear 
medium and (ii) conventional nuclear-structure effects, such as binding-energy corrections 
and Fermi motion. The first relativistic description of neutrino-nucleus scattering address- 
ing these complications was presented in Ref. where the authors employed a relativistic 
Fermi- gas model for the target nucleus. Nuclear-structure effects were further investigated 
in Ref. [15| by considering, in addition to a relativistic Fermi-gas model, a description of the 
bound nucleon in terms of harmonic oscillator wavefunctions. In general it was found that 
nuclear-structure effects and final-state interactions can have a substantial effect on the 
individuai cross sections. However, as originally suggested in Ref. 0], the ratio of proton 
to neutron yields is largely insensitive to these effects 0,0. 

In this work we also present a relativistic description of neutrino-nucleus scattering, but 
employing bound-state wavefunctions obtained from an accurately calibrated mean-field 
model [19j . Further, the impulse approximation is assumed, that is, the neutrino-nucleon 
interaction is assumed unchanged in the nuclear medium. Finally, as the ultimate aim of this 
project is to compute ratios of cross sections, distortion effects on the ejectile nucleon will 
be neglected. As we will show later, this leads to a great simplification in the calculation of 
ali relevant quantities. The paper has been organized as follows. The formalism is presented 
in Sec. HH followed by results and conclusions in Sec. Hill and Sec. IIV[ respectively. 



II. FORMALISM 

In this section the formalism for the relativistic description of (neutral-current) neutrino- 
nucleus scattering will be presented. In particular, it will be shown in Sec. Ili Al that the 
cross section can be written as a contraction between leptonic and hadronic tensors. In turn, 
by relying exclusively on fundamental principles, the hadronic tensor will be decomposed in 
terms of a set of invariant structure functions (see Sec III B|) . Thus, the formalism is model- 
independent. Yet to determine the structure function, and ultimately the cross section, one 
must rely on a model. This will be discussed in Sec. Ili CI 



A. Cross section in terms of leptonic and hadronic tensors 

Due to the short-range nature of the weak interaction, the one-boson exchange approx- 
imation provides an excellent description of neutrino-nucleus scattering. This results in a 
cross section that cleanly separates (or factorizes) into leptonic and hadronic components. 
The kinematics of the process is depicted in Fig. ^ Here the initial and final neutrino 
four-momenta are denoted by k and k', respectively. Further, the reaction proceeds via 
the exchange of a virtual Z°-boson with four-momentum q. The target and residuai nu- 
cleus have four-momenta denoted by P and P', respectively. Finally, the ejectile proton has 
four-momentum p' and spin component s'. 

The differential cross section can now be defined in terms of these kinematic variables 
and the transition matrix element Ai as follows: 

(2tt) 4 5{k + P - k' - p' - P' ) d 3 k' d 3 p' d 3 P' 2 

| Vl - v 2 | (2^ (2^ (2^p 1 1 ' ( ) 

where v x — v 2 denotes the initial relative velocity. The transition-matrix element contains 
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FIG. 1: Lowest-order Feynman diagram for the knockout of a bound nucleon via neutral-current 
neutrino-nucleus scattering 

ali the dynamical information about the reaction and is given, in the conventions of Ref . , 
by 



4M 
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Note that in Eq. (J2J) the initial and final nuclear states are denoted by ^i(P) and tyf(P'), 
respectively. Further, g is the weak coupling Constant and J^(q) is the weak nuclear current 
operator, which contains both vector and axial-vector components. Finally, as only low 
momentum transfer (\q 2 \ <C M|) scattering will be considered, the following replacement is 
valid: 



D»„{q) 



9/iu + q^v/ M z 



q 2 -M 2 z 

This allows the transition matrix element to be written as 
Gp 



9jw_ 
M z 
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Note that Gf is the Fermi Constant for muon decay which is given by 

G p 



A^2M 2 



1.166 x 10 GeV~ 



(3) 
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In Eq. (J2J) and throughout this work plane-wave Dirac spinors are defined as follows: 

(E k = Vk 2 + M 2 



W(k, s) 



E k + M 



2Eh 



1 

cr • k 



Xs 



Ek + M 



(6) 



Note that the above defìnition corresponds to the normalization 



(7) 



This (non-covariant) normalization is motivated by the standard choice adopted for bound- 
state spinors (see Sec. IH C| and Ref. |20j) which is given by 



j Ul{v)U a {r) dx=\. 



(8) 



For (massless) neutrinos in the laboratory frame the initial fìux factor in Eq. (JTJ is equal to 
one. Substitution of the above relations into Eq. (JTJ leads to the following expression for 
the differential cross section: 



r 2 r/ 3 V W 3 tV ^7 3 P' 

da = ^f(2^5(k + P — k' —p' — P')44 „ 44 LrW" r 



(2vr) 3 (2vr) 3 (2tt) 
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where the leptonic tensor is given by 



F(k') (7" -7V M k ) F(k') (7" -7 i; 7>(k) 



while the hadronic tensor by 



(p',s'; Vf(P') >(g) %(P) (p',s'; * f (P') r(q) %(P) 



(9) 



(10) 



(11) 



The integrai over P' may be performed using the spatial part of the Dirac delta function. 
This fixes the three-momentum of the residuai nucleus to be 



P' = k - k' - p + P > q - p' 



lai) 



(12) 



where q = k— k' is the three-momentum transfer to the nucleus. The differential cross section 
can now be written as 



da 



G 



1 d 3 k'd 3 p'6(E k +M A -E k/ -E pl -E pl )L u W^ . 



K2vr) ; 



(13) 



B. Differential cross section in terms of nuclear structure functions 



In Sec. Ili Al it has been shown that the differential cross for nucleon knockout in neutrino- 
nucleus scattering involves the contraction between the leptonic tensor [Eq. ((TUJ)] and the 
hadronic tensor W^ u [Eq. (jll|) ]. In this section both of these quantities will be calculated 
in a model-independent way by introducing a suitable set of nuclear structure functions. 

Starting from Eq. ()1U|) it follows that the leptonic tensor may be written as 



Tr 
2 

kk' 



( 7 M_yy) ^(k)F(k) j (/- 7 V) (W) 

k»k w + k'»k v - g^k ■ k' + ihe^ al3 k a k'p 



(14) 
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where we made use of Eq. ((HJ) to write 



(*=|k|). (15) 



Note that in the above expressions /i = — 1 and h — +1 refer to left-handed neutrinos and 
right-handed antineutrinos, respectively. For later convenience, the leptonic tensor can now 
be separated into a symmetric and an antisymmetric part. That is, 



£^ = £^+£^, (16) 



where 



1% = — (Wk' v + k'»k u - g^k ■ k!) , (17a) 
kk' 

C = j^ihe^Kk'p . (17b) 

Note that ali remnants of the neutrino helicity resides in the terni containing the antisym- 
metric Levi-Civita tensor. It is only this antisymmetric component of the leptonic tensor 
that is sensitive to the difference between an incident neutrino or antineutrino beam. It 
then follows from Eq. ()17|) . as the weak neutral currents is conserved for massless neutrinos, 
that 

% l*> = £^q u = , (18) 

where g M = (lj, q) = (k^ — k ,fl ) is the four-momentum transfer to the nucleus. 

The hadronic tensor is an extremely complicated object as in principle exact many- 
body wave functions and operators must be used. Yet it follows from Eq. (|TT|) that for 
unpolarized nucleon emission the hadronic tensor is only a function of three independent 
four-momenta: P M and as the four-momentum of the recoiling nucleus P' M is fixed 
by four-momentum conservation. The hadronic tensor can therefore be parametrized in 
terms of a basis constructed from the following five tensors: ^q tl ,P^,p" tJ, ,g tJ,l ',e lll ' a ^^. This 
is analogous to the case of electron scattering but now we are no longer allowed to invoke 
either parity invariance or current conservation, as the weak interaction violates parity and 
the axial-vector current is not conserved. We start by separating the hadronic tensor into 
symmetric and antisymmetric components. That is, 

W*» = Wg v + . (19) 

Using the above-mentioned basis the individuai components may be written as follows: 

Wg v = Wxg^ + W 2 q^q u + W 3 P> M P U + W A p ,fM p' u 

+ W 5 (q^P v +P^q v ) + W 6 (q^p' u + p^q v ) + W 7 (P^p' u + p'^P v ) , (20a) 

= Wz{q»P v -P»q u ) + W 9 {q»p' u - p'»q u ) + W w {P»p' v - p'»P v ) 

+ W u e^q a Pp + W 12 e^ aP q a p' p + W 13 e^ P a p' p . (20b) 

Note that ali structure functions are functions of the four Lorentz-invariant quantities, q^q^ = 
— Q 2 , q- P, q-p' -, and P-p'. The hadronic tensor for the (u, v'p) [or (u, v'n)] reaction contains 
thirteen independent structure functions. Contrast this, for example, to the hadronic tensor 
for (e, e'p) which is fully written in terms of only five structure functions 
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We now proceed to evaluate the contraction of the leptonic tensor with the hadronic 
tensor. First, we introduce the following definition: 



T{k, k'; P,p') ee (j^J £, V W^ = T s {k, k'; P,p') + T A {k, k'; P,p') 



(21) 



where an d J~A are defined in terms of the symmetric and antisymmetric components of 
the leptonic and hadronic tensors, respectively. That is, 



F s (k,k';P,p') 



F A (h, k';P,p') 



4 

kk' 



9 s w ^ 
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(k.P)(k'-P) + ^ (k-k') 



+ w 4 



+ w- 



M 2 

(k ■ P')(k' ■ p') + -f(k ■ k') 



(k ■ P)(k' ■ p') + (k ■ p')(k' -P)-(k- k')(P ■ p') 



4 

W 



-i 



9 A W iw 



= ih \W 10 e^kXP a p'p + Wu(k -k')(k-P + k'- P) 
+ W l2 (k- k')(k-p+k'-p 



+ W 



13 



{k-p'){k' -P) - {k-P){k' -p') 



(22a) 



(22b) 



Note that as a result of current conservation, only eight of the originai thirteen structure 
functions survived the contraction [see Eq. (j!8|) ]. Further, an interesting difference between 
neutrino and electron scattering can be seen from Eq. (j22b|) . While for electron scatter- 
ing one must prepare a polarized beam to sample the antisymmetric part of the hadronic 
tensor j2l|, for neutrinos the polarization happens by default. Eqs. (|22a|) . and ()22b|) 

comprise the principal results for this section. It is the most general structure possible for 
nucleon knockout in neutral-current neutrino-nucleus scattering. It shows that the differen- 
tial cross section is completely determined by a set of eight structure functions multiplied 
by kinematical factors. 



C. Model-dependent calculation of the cross section 



In the previous section it was shown that the differential cross section is completely 
determined by a set of eight structure functions. These structure functions parametrize our 
ignorance about strong-interaction physics. In principle, these structure functions could be 
measured through a "super" Rosenbluth separation. In practice, however, this is beyond 
realistic expectations. Thus, it is not possible to proceed further without an explicit model 
of the hadronic vertex. 



7 



First, we focus ori some "kinematical" approximations that are macie in order to simplify 
the argument of the energy conserving delta function in Eq. ()13p. In the laboratory frame 
the total energy of the residuai nucleus is given by 

E P , = ^P^ + = ^(q - p') 2 + M\_ x « M A _ X , (23) 

where the last approximation follows in the limit of no recoil corrections. Further, it is 
assumed that the energy transfer to the nucleus satisfies 

u = k — k' = E p i — Eb , (24) 

where Eb is the energy of the struck nucleon. This implies that 

E P , = M A -E B , (25) 

and justifies the following replacement in Eq. (JTSjl : 

5(E k + M A -E k ,-E p ,-E P ,) — > 5(k~k'-E p/ + E B ) . (26) 

Note that we have defined E^— |k| = k and Ey = |k'| = k! . 

Z A 




FIG. 2: Coordinate axes used to define the angles a, 9 and 4>. 

Consider now the coordinate system shown in Fig.|21 where the incoming neutrino defines 
the z-axis. The outgoing nucleon is detected at a scattering angle a relative to the z-axis, 
while the outgoing neutrino, assumed undetected, has polar and azimuthal angles 9 and 0, 
respectively. The following relations are therefore valid: 

d 3 k' = k' 2 sm 9dk'd9d0 and d 3 p' = 2>!rp'E p/ dE p ,d(cosa) . (27) 

As the outgoing neutrino will remain undetected, one must integrate over k', 9 and <fi. Using 
Eqs. (j2IIl) and (J2Tj) the following expression for the differential cross section [Eq. (fT3j)] is 
obtained: 

IeM^) = &) (") l'^C"W;P,v') , (28) 



s 




FIG. 3: Graphical representation of the approximation employed at the hadronic vertex to obtain 
a tractable form for the matrix element given in Eq. (|29|) . 



where the energy conserving delta function constrains the energy of the outgoing neutrino 
to \k'\ = \k\+E B -E p ,. 

What remains now is to provide an explicit form for the eight independent structure 
functions introduced in the previous section. It follows from Eq. (0J) that the dynamical 
information on the hadronic vertex is contained in the following matrix element (and its 
complex conjugate): 

J»={p',s';* f (P')\>(q)\%(P)) . (29) 

To obtain a tractable form for this extremely complicated object we rely on the approxima- 
tions depicted in Fig.El that we now address in detail. First, it is assumed that the Z°-boson 
couples to a single bound nucleon. This neglects two- and many-body components of the 
current operator. Second, it is assumed that the detected nucleon is the one to which the 
virtual boson couples to. This neglects two- and many-body rescattering processes. Finally, 
we neglect final-state interactions (distortions) of the ejected nucleon. While many similar 



treatments incorporate distortion effects on the ejectile (see for example Ref. [21| in the 
case of electron scattering) concentratine on ratios of cross sections makes the formalism 
largely insensitive to distortion effects [12]. This will render the hadronic tensor analytic. 
Incorporating the above simplifications yields a hadronic matrix element of the following 
simple form: 

J" = W(p / ,s / )J' i (?)Z4(p m ) , (30) 

where a represents a collection of quantum numbers and p m = p' — q is the missing momen- 
tum. Note that the bound-state wavefunction is given by 

u*(p) = zwp) = -h)' ( f J$$. P ) ì y™<p) • ( 31 ) 



Here E is the bound-state energy, k the (generalized) angular momentum, and m the spin 
projection. Further, QekÌJ?) and Jek(p) are Fourier transforms of the upper and lower com- 
ponents of the bound-state wavefunction, respectively |Ì3 |. 
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The impulse approximation is now invoked by assuming that the weak neutral current 
for a nucleoli in the nuclear medium retains its free-space form. That is, 



J„ = J, 



NC 



9NC 
"V5 



F 1 (Q 2 ) lfl + iF 2 (Q 2 )a 



fJtU 



2M 



(32) 



Here Fi and F 2 are Dirac and Pauli vector form factors, respectively and Ga is the axial 
form factor. The additional induced pseudoscalar form factor (proportional to q^) does 
not contribute to (massless) neutrino scattering and will be neglected henceforth. The two 
vector form factors may be decomposed in terms of the usuai electromagnetic Dirac and 
Pauli form factors, plus a yet undetermined strange-quark contribution. Similarly, the axial 
vector form factor consists of a purely isovector contribution, that may be determined from 
Gamow-Teller /3-decay rates, and a purely isoscalar strange-quark contribution. It is the aim 
of this contribution to explore the sensitivity of neutral-current neutrino-nucleus scattering 
to the strange-quark contribution to the axial form factor. A detailed discussion of the weak 
neutral current [Eq. (|3*2*)l ] is given in Appendix 1X1 

It is now possible, using Eq. ()30p. to explicitly calculate the hadronic tensor defined in 
Eq. (fTT|). Assuming that the spin of the outgoing nucleon is not detected, the hadronic 
tensor may be written as 



(33) 



where J v = 7°( J")^ . The normalization used in Eq. (J7J) implies that the standard (on-shell) 
Feynman propagator is given by 



= £^[W(p^0>(<Z)W a , m (p m ) U(p', S ')r(q)U a , m (p r , 

s' m 

= Tr \j^q)S a {p m )r(q)S{p') 



5(p') = ^W(p', S ')W(p / ,, / ) + 



2E P > 



(34) 



Moreover, it has been shown in Ref. 2^1 that the following simple identity is valid even in 
the case of a bound-state spinor: 



S a (p) = Ys^Up) U a , m {p) = (2j + 1)(^ + M a ) 



(35) 



Note that in the above equation for the "bound-state" propagator mass-, energy-, and 
momentum-like quantities have been introduced. These are given by 



Pa 



'gì(p) + fl(p) 

29a(p)fa(p)p 



(36a) 
(36b) 
(36c) 



and satisfy the "on-shell relation" 



pì = eI-pI = m ( 



(37) 
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It now follows from Eqs. and f)35[) that the hadronic tensor may be written in the 

following simple forni: 



2J + 1 
2E rt i 



Tr 



>( g )(^ + M a )J^)(/ + M) 



(38) 



The fact that the hadronic tensor may be expressed as a trace over Dirac matrices, despite the 
presence of the bound-state wave function, greatly simplifies the calculation. It is important 
to note, however, that this enormous simplification would have been lost if distortion effects 
would have been incorporated in the propagation of the emitted nucleon. The emphasis on 
computing ratios of cross sections is the main justification behind this simplification. We 
now obtain 



2J + 1 
2E P , 



Wi<T + W 2 (p'»p u a + + W 3 (« + p£q u ) + W 4 W + pV) + 



W^q v + W 6 e^ x p a ^p' x + W 7 e^ x p a ^q x + W 8 e^ x p' a q x 



(39) 



Ali components of the hadronic tensor are given in terms of relatively simple expressions 
that have been collected in Appendix [Bj These model-dependent structure functions may 
be related to the model-independent ones [see Eqs. (J2Ua|) and (|20b|) ] as has been shown in 
Appendix [HJ This concludes the formalism for neutral-current neutrino-nucleus scattering. 
The explicit forms of the structure functions given in the appendix can now be used to 
evaluate the differential cross section defined in Eq. 



III. RESULTS 



In Ref. H3 the strange-quark content of the nucleon was studied via neutrino-nucleus 



scattering by using a relativistic Fermi gas model of the target nucleus. This amounts 
to averaging the free neutrino-nucleon cross section over a sharp momentum distribution 
for the struck nucleon. In this work we improve on the above description by employing 
bound-nucleon wavefunctions obtained from a relativistic mean-field approximation to the 
accurately calibrated NL3 model of Ref. 0|. To quantify the impact of this improvement we 
display in Figs. lUtoElthe doublé differential cross section d 2 a / dE p >d(cos a) as a function of 
the ejectile nucleon kinetic energy T p and its scattering angle a in the laboratory frame. The 
incident neutrino energy has been fixed in these plots at 150, 500 and 1000 MeV, respectively. 
For illustration purposes — and only for these three graphs — the strange-quark contribution 
to the axial-vector form factor (g\) has been neglected and only proton knockout from the 
lp3/2 orbita] of 12 C is considered. 

For the lowest energy neutrinos the cross section displays a single well developed peak 
at low T p that monotonically decreases with increasing scattering angle a. For E v = 500 
MeV, our Fig. |5] may be compared directly to Fig. 1 of Horowitz and collaborators [1J| . In 
particular, for a scattering angle of a = 20°, our calculation (third curve along the a direction) 
also exhibits the characteristic double-humped structure. For larger values of a the peaks 
merge into one and the cross section develops a shape similar to that of E u = 150 MeV. Yet, 
an important difference between the two sets of calculations is that our cross section does 
not develop the sharp features displayed at small angles in Ref. This is due to the more 
realistic momentum distributions used in our calculations. 
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d 2 o/dE d(cos a) (fm 2 /MeV) 




FIG. 4: Doublé differential cross section d 2 a/dEd(cos a) as a function of the outgoing proton 
laboratory kinetic energy T p and laboratory scattering angle a. The calculation shown is for 
proton knockout from the lp 3 / 2 orbitai of 12 C at an incident neutrino energy of 150 MeV. 



Next, we produce angle-integrated cross sections as a function of T p . Figs.0 and |H1 display 
cross sections for the knockout of protons and neutrons, respectively. The long-dash-short- 
dashed line represents knockout from the ls 1 / 2 orbitai while the dashed line from the lp 3 / 2 
orbitai of 12 C; the solid line displays their sum. As before, calculations are shown for incident 
neutrino energies of 150, 500, and 1000 MeV, respecti vely . The cross sections at 150 MeV 
correspond to those shown in Figs. 4 and 5 of Ref. jl4|. It has been shown in Ref. 
that binding-energy corrections (at 150 MeV) reduce the cross sections relative to their free 
Fermi-gas values by about 40%. Our cross sections (already summed over both occupied 
orbitals) are reduced even further. Note that while an average binding energy of 27 MeV has 
been used in Ref. |l4j . our calculations include binding energies computed exactly within a 
mean-field approach. At the higher energies of 500 and 1000 MeV the cross sections display 
the same general trend as the one for 150 MeV, namely, a peak at a low value of T p and a 
"smooth" falloff with increasing T p . 

In Figs. 01 and we examine the impact of a strange-quark contribution ( q s a) to the 
axial form factor on the cross section. Following the discussions in Refs. 0, 0, 14]| a 
value of g% = —0.19 is adopted henceforth. Further, in ali cases presented here strange- 
quark contributions to the weak vector form factors are ignored. The solid and dashed 
lines correspond to a zero and a non-zero value of g s A , respectively. In both figures we have 
summed over the ls 1//2 and the lp 3 ^ 2 orbitals of 12 C. Because of the dominance of the axial- 
vector form factor (see discussion below) a non-zero value of g s A increases the cross section 
for proton knockout, whereas for neutron knockout the cross section is decreased. These 
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FIG. 5: Doublé differential cross section d 2 a/dEd(cos a) as a function of the outgoing proton 
laboratory kinetic energy T p and laboratory scattering angle a. The calculation shown is for 
proton knockout from the lp 3 / 2 orbitai of 12 C at an incident neutrino energy of 500 MeV. 

findiners are consistent with those of Refs. 

Next we investigate the role of the various single-nucleon forni factors in the calculation 
of the differential cross section. For this case we restrict ourselves to proton knockout from 
the lp 3//2 orbitai of 12 C. As before, we consider incident neutrino energies of 150, 500, and 
1000 MeV. The results are shown in Fig. E] in this figure the solid line corresponds to 
the case (G A ^0,Fi^ 0, and F 2 ^ 0), the dashed line to (G A ^ 0, F x = 0, and F 2 = 0), 
the long-dash-short-dashed line to (G A ^ 0, Fi = 0, and F 2 ^ 0), and the dash-dot line to 
(G A = 0, Fi 7^ 0, and F 2 ^ 0). For ali energies we observe the dominant role played by the 
axial-vector forni factor Ga- Indeed, for 150 MeV there is virtually no distinction between 
the calculation using non-zero values for ali forni factors (solid line) and the one where only 
the axial-vector forni factor is included (dashed line). This is due to the smallness of the 
weak vector charge of the proton (at ali values of Q 2 ) and the low-momentum transfer of the 
reaction, which makes the contribution from the Pauli form factor small. The dominance 
of the axial-vector form factor is further illustrated by the fact that when it is set to zero, 
the cross section becomes vanishingly small (dash-dot line). This behavior is important as 
it increases the sensitivity of the reaction to the strange-quark contribution to the the axial 
form factor. Indeed, for Fi = F 2 = 0, the differential cross section becomes proportional to 
the square of the axial-vector form factor which, at Q 2 = 0, is given by (see Appendix |A"|) 

G 2 A (Q 2 = 0) = (g A - g s A f = (g A + {g s A f - 2g A g s A ) • (40) 
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FIG. 6: Doublé differential cross section d 2 a/dEd(cos a) as a function of the outgoing proton 
laboratory kinetic energy T p and laboratory scattering angle a. The calculation shown is for 
proton knockout from the lp 3//2 orbitai of 12 C at an incident neutrino energy of 1000 MeV. 



The sensitivity to the strange forni factor comes about through the interference terni. 

An important problem encountered in the previous (and most of the earlier) analysis is the 
dependence of the cross section on distortion effects, that is, on the final-state interactions 
between the ejectile nucleon and the residuai nucleus [l5l . ufy. In an effort to circumvent 
this problem, the authors of Ref. [TtJ proposed to measure the ratio of quasielastic proton 
yield to quasielastic neutron yield, rather than the absolute cross section. Focusing on the 
ratio of cross sections proves to be advantageous for a number of reasons. For example, 
the calculation of the angle-integrated cross section in Ref. Q is particularly sensitive to 
binding-energy corrections: both proton and neutron knockout cross sections are quenched 
by almost 40% relative to the free Fermi-gas estimate. Unfortunately, a strange-quark con- 
tribution to the axial-vector form factor of the neutron also reduces the cross section by 
nearly 40%. Hence, it might be dimcult to separate nuclear-binding effects from a genuine 
strange-quark contribution. (We reiterate that, while it remains advantageous to reduce the 
sensitivity of the reaction to nuclear-structure effects, the merit of our calculation is that it 
incorporates realistic binding energies and momentum distributions.) Further, while distor- 
tion effects modify the cross section, they often do so without a significant redistribution of 
strength. Thus the ratio of cross sections, rather than the cross sections themselves, should 
be less sensitive to distortion effects. Indeed, in the model of Ref. 3| it was shown that the 
ratio of cross sections was insensitive to distortion effects. Hence, we conclude this section 
by displaying in Fig. ^Jthe ratio of cross sections for protons over that for neutrons. While 
for the lowest value of the neutrino energy (150 MeV) the ratio remains fairly Constant, a 
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FIG. 7: Differential cross section da/dE as a function of the outgoing proton laboratory kinetic 
energy T p . The dashed and long-dashed-short-dashed lines are for proton knockout from the lp 3//2 
and ls 1 / 2 orbitals of 12 C, respectively; the solid line represents their sum. The incident neutrino 
energy is taken to be E k = 150, 500 and 1000 MeV. 



signifìcant dependence on the outgoing nucleoli kinetic energy T p is observed for the other 
two cases. How sensitive this dependence is to the high-momentum components in the nu- 
clear wavefunction (induced, for example, by short-range correlations) remains an important 
open question. 



IV. SUMMARY AND CONCLUSIONS 



Neutral-current neutrino-nucleus cross sections have been computed in a relativistic piane 
wave impulse approximation. The aim of this contribution was to examine the sensitivity of 
the reaction to the strange-quark contribution to the axial-vector forni factor of the nucleon. 
This was done, in part, in response to the Fermilab Intense Neutrino Scattering Experiment 
(FINeSE) initiative. In this work we have followed closely the seminai contributions made to 
this subject by various groups [14. 

E3)E3|- Yet, we have improved significantly on them by 
incorporating nuclear structure effects through an accurately-calibrated relativistic mean- 
field model. Thus, accurate binding energies and nucleon momentum distributions were 
employed. 

Our results indicate signifìcant quantitative, although minor qualitative, changes with 
respect to the relativistic Fermi-gas model of Horowitz and collaborators First, double- 
differential cross sections displaying sharp features due to a discontinuous (Fermi gas) mo- 
mentum distribution get soften by our choice of mean-field wave functions. However, most 
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FIG. 8: Differential cross section da/dE as a function of the outgoing neutron laboratory kinetic 
energy T p . The dashed and long-dashed-short-dashed lines are for neutron knockout from the lp 3 / 2 
and ls 1 / 2 orbitals of 12 C, respectively; the solid line represents their sum. The incident neutrino 
energy is taken to be E k = 150, 500 and 1000 MeV. 



of the sharp features of the Fermi-gas cross section disappear as soon as one integrates over 
the scattering angle of the ejected nucleon. At this point the shape of the cross section 
becomes largely insensitive to the choice of momentum distribution. Not so its magnitude. 
Low-energy cross sections with binding-energy corrections included in an average way were 
shown to be reduced by 40% relative to the corresponding free Fermi-gas estimates • 



The cross sections reported here, computed with binding energies obtained from a relativis- 
tic mean-field model, yield even smaller cross sections. 

In order to reduce the sensitivity of the strange-quark content of the nucleon to nuclear- 
structure effects, while at the same time eliminating the need for an absolute normalization 
of the cross section, the authors of Ref. [ItJ have proposed to use the ratio of proton to 
neutron yields. Indeed, these authors demonstrated that the ratio of cross sections is largely 
insensitive to distortion effects [l6[. Further, binding-energy corrections, which were respon- 
sible for the large (40%) reduction of the absolute cross section, become barely visible in 
the cross-section ratio 14]. In our case the insensitivity to final-state interactions entailed 



an enormous simplification: by introducing the notion of a bound-nucleon propagator we 
have exploited Feynman's trace techniques to develop closed-form, analytic expressions for 
the cross section. 

In summary, the sensitivity of neutral-current neutrino-nucleus scattering to the strange- 
quark content of the axial-vector form factor of the nucleon was examined. A model- 
independent formalism based on a set of eight nuclear structure functions was developed. 
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FIG. 9: Effect of a strange-quark contribution (of g s A = —0.19) to the axial-vector form factor on the 
differential cross section da/dE as a function of the laboratory kinetic energy T p of the outgoing 
proton. The solid and dashed lines correspond to a zero and a non-zero value of g A , respectively. 
In this figure we have summed over the ls 1 / 2 and Lp 3 / 2 orbitals of 12 C. 



On account of both, the notion of a bound-state nucleon propagator and the insensitiv- 
ity of the ratio of proton-to-neutron yields to distortion effects, we computed ali nuclear 
structure functions in closed form. Adopting a value for strange-quark contribution to the 
axial-vector form factor of the nucleon of g s A = —0.19, led to a significant enhancement in 
the proton-to-neutron yields relative to one in which the strange-quark contribution was 
neglected. 
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FIG. 10: Effect of a strange-quark contribution (of g s A = —0.19) to the axial-vector forni factor 
on the differential cross section da/dE as a function of the laboratory kinetic energy T p of the 
outgoing neutron. The solici and dashed lines correspond to a zero and a non-zero value of g s A , 
respectively. In this figure we have summed over the ls 1 / 2 and lp 3 / 2 orbitals of 12 C. 

APPENDIX A: HADRONIC WEAK-NEUTRAL CURRENT 

Neutral current vector and axial-vector matrix elements between on-shell nucleon states 
can be parametrized on general grounds in terms of four forni factors in the following forni: 



<iV(p', S ')|j; c (0)|iV(p, s ))= M (p', s ' 



F 1 {Q 2 ) lll + iF 2 {Q 2 )a i 



ni/ 



2M 



u(p,s), (Ala) 



(N(p>,s>)\J^(0)\N(p,s))=u(p',s>) [G A (Q' 2 h, + G P (Q 2 )^\ 76 «(p, a ) . (Alb) 

The vector current may be decomposed in terms of isoscalar and isovector electromagnetic 
contributions plus an explicit strange-quark contribution, which is assumed isoscalar. Simi- 
larly, the axial-vector current may be related to the isovector current measured in neutron 
beta decay plus an isoscalar strange-quark contribution. That is, 

J* c = (2 - 4 sin 2 9 W ) JJ M (T = 1) - 4 sin 2 9 W J™(T = 0) - s^s , (A2a) 
J„ N 5 C = -2À,(T=1) + s W s . (A2b) 
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FIG. 11: Effect of the single- nucleon forni factors on the differential cross section da/dE as a 
function of the laboratory kinetic energy T p of the outgoing proton. The calculation shown is for 
proton knockout from only the lp 3 / 2 orbitai of 12 C at incident neutrino energies of 150, 500 and 
1000 MeV. Explanation for the various lines is given in the text. 



This decomposition enables one to express the two vector and the one axial-vector forni 
factors in the following forni: 

%{Q 2 ) = (1 - 4sin 2 fl w )i 7 ? ) (Q 2 ) - F*> n \Q 2 ) - F^ S \Q 2 ) ; forprotons , (A3a) 

F^Q 2 ) = (1 - 4sin 2 fl w )if l) (Q 2 ) - if^Q 2 ) - if } (Q 2 ) ; forneutrons , (A3b) 

and 

G A (Q 2 ) = +Gf(Q 2 ) - G { °\Q 2 ) ; forprotons , (A4a) 

G A {Q 2 ) = -Gf{Q 2 ) - G ( a\Q 2 ) ; forneutrons . (A4b) 

Note that because (massless) neutrino scattering is insensitive to the induced pseudoscalar 
forni factor Gp, it has been ignored throughout this work. In what follows, standard pa- 
rameterizations of the electromagnetic and axial-vector nucleon forni factors are employed. 
In particular, we follow the conventions adopted in Ref. 0. These are given by 

F} P \Q 2 ) = ( 1+ ; ( ^ Ap) ) G V D (Q 2 ) , F^(Q 2 ) = G V D (Q 2 ) , (A5a) 

F[ n \Q 2 ) = ( AnT 1 ( ^; ?7) ) GUQ 2 ) , Ft\Q 2 ) = ( An( 1 1 ^ 7 " ?) ) ' (A5b) 

Gf{Q 2 ) = 9a GUQ 2 ) , G^(Q 2 ) = g s A G A D {Q 2 ) , (A5c) 
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FIG. 12: Ratio of cross sections for protons over that for neutrons as a function of the laboratory 
kinetic energy T p of the outgoing nucleon. Cross sections include contributions from the ls 1 / 2 and 

l p 3/2 orbitals of 12 C _ 



where a dipole forni factor of the following forni is assumed: 

G V D (Q 2 ) = (1 + Q 2 /M*)- 2 = (1 + A.97ry 2 (A6a) 

GUQ 2 ) = (1 + Q 2 /M 2 A )- 2 = (1 + 3.31r)- 2 (A6b) 

rj = (1 + 5.6 t)- 1 t = Q 2 /{AM 2 ) . (A6c) 

Finally, for reference we display the value of the various nucleon forni factors at Q 2 = 

F 1 (p) (0) = l, F 1 (n) (0) = 0, F 1 (s) (0) = 0, (A7a) 

^ 2 (P) (0) = A P = +1.79, F 2 W (0) = A n = -1.91 , F 2 (s) (0) = /x a , (A7b) 

Gf(0) = ^ = +1.26, G { °\0) = g s A . (A7c) 



APPENDIX B: HADRONIC TENSOR IN A RPWIA 



The hadronic tensor computed in a relativistic plane-wave impulse approximation has 
been written in Eq. (J39)) in terms of eight structure functions. These are given by the 
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following simple forms: 
Wi 



AF 2 {MM a -p a . p ') + -^ \{MM a + p a ■ p')q 2 - 2{p a ■ q)(p' ■ q) 



- AG 2 A (MM a +p a -p') + 



4iqF 2 
M 



(p' ■ q)M a - (p a ■ q)M 



W 2 

w 3 
w 4 

w 5 

w 6 
w 7 



A[F( + G A ) l- 2 

^2 fp'-v 



2[jr_ 

M 2 



M 2 



2 fPa-q 



M 2 



+ 2F 1 F 2 , 



F 2 



M 2 
ZìFxGa 
AiF 2 G a 



(MM a +p a -p') 



4i 



M 



F 2 Ga 



(Bla) 
(Blb) 

(Blc) 

(Bld) 

(Ble) 

(Blf) 
(Big) 

(Bili) 



Our final task consists in relating the structure functions in the general model- 
independent expansion of W^ v given in Eqs. ()20a|) and (|20b|) . to the model-dependent ones 
given above. This is done by noting that in the laboratory frame one can write 



p^ a = aP» + b{p'^ - g") 
where the a and b coefficients are defined as follows: 



M A 



E a 



-En 



Pr. 



IPc 



Pr. 



(B2) 



(B3) 



Substitution of Eq. ()B3|) into Eq. (|3*9J) allows us to identify the contribution of Wi to each 
of the structure functions. These are given by 



W 4 
W 7 

w 13 



2j + l 
2E p > 

2J + 1 
2E P , 

2J + 1 
2E p , 

2J + 1 
2E P > 

2j + l 
2E p / 

2J + 1 
2E p i 





(B4a) 


(2bW 2 ) , 


(B4b) 


(aW 2 ) , 


(B4c) 




(B4d) 


(-bW 6 -bW 7 -W 8 } , 


(B4e) 




(B4f) 
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Note that the model predicts W3 = W\o = 0. Moreover, as was mentioned previously, 
(massless) neutrino scattering is insensitive to the following five structure functions: W2, 
W 5 , W 6 , W 8 , and W 9 . 
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